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Abstract 

We study the long time behavior of a Brownian particle moving in 
an anomalously diffusing field, the evolution of which depends on the 
particle position. We prove that the process describing the asymptotic 
behaviour of the Brownian particle has bounded (in time) variance when 
the particle interacts with a subdiffusive field; when the interaction is 
with a superdiffusive field the variance of the Umiting process grows in 
time as t^^~^ , f/2 < 7 < f. Two different kinds of superdiffusing (ran- 
dom) environments are considered: one is described through the use of 
the fractional Laplacian; the other via the Riemann-Liouville fractional 
integral. The subdiffusive field is modeled through the Riemann-Liouville 
fractional derivative. 

Keywords: anomalous difFusion, Riemann-Liouville fractional derivative (inte- 
gral), fractional Laplacian, continuous time random walk. Levy flight, scaling 
limit, interface fluctuations. 



1 Introduction 

In [1], L. Bertini et al considered the following system of Ito-SDEs, describing 
the evolution of a one-dimensional interface: 

dX{t) = \dw(t) + a{ipx(t).h{i))dt 

(1) 

dh{t) = \/^h{t)dt ~ ipx{t)dX{t), 

with initial conditions X{Q) = h{0) = 0. In the above system w{t) is a one di- 
mensional Brownian motion (BM) on the filtered probability space (fi, Tt,P) 
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{E is going to denote expectation with respect to P) and (•,•) is the scalar 
product of L^(M,dx). More precisely, in [1] the authors consider a system 
thermally isolated from the exterior, in a state in which two phases coexist. 
Under the assumption of planar symmetry for the system, the interface po- 
sition is represented by the point X{t) G C(M+) separating the two phases. 
In equation ([T])i the interface displacements are described as the sum of two 
contributions: the first is a Brownian fluctuation, related to the macroscopic 
fluctuations of the system, the second is the interaction with a diffusive field, 
h{t) = h{t,x) g C(]R+;C(R)). Also, 

(ifxit) , Ht)) = ^ dxifiix - X{t))h{t, x), 

where ^p{x) is a probability density in the Schwartz class (regions of the field 
far from the interface do not significantly affect the interface evolution) and 
Vx(t) = ^p{x- X{t)). 

On the other hand, equation ^2 describes the field variation as the sum of a 
diffusive term plus a "feedback term" taking into account the latent heat effect. 
The parameters A > and a > determine the intensity of the Brownian noise 
and of the coupling with the field, respectively. In [I] the authors study a scaling 
limit of X{t) as A — >■ under the hypothesis a = A of weak coupling. 
Notice that the system ([TJ can also be interpreted as describing a Brownian 
motion weakly coupled with a (diffusive) random environment, the evolution of 
which depends on the position of the Brownian motion itself. For further details 
about the model we refer to [IJ and |2j. 

Let ^(t) be the solution of the following integral equation 

^{t)^h{t)~ fdspt-smis), (2) 

J o 

where b{t) is the scaled BM b{t) — Xw{tX^^) and pt{x) = p{t,x) is the density 
of a centered Gaussian with variance t. In [Ij the following asymptotics ([3]) 
and Q are obtained: upon rescaling the interface position, i.e. considering the 
process Xx{t) = X{tX-'^), we have that \/N e [1,(x))3t = t{N) > s.t. 

limE sup \Xx{t) - ^{t)\^ ^0. (3) 

-^^0 t<T|logA 

As noticed in [I , this implies that X\ converges weakly to ^ as A — > in C(K+) 
endowed with the topology of uniform convergence on compacts. Furthermore, 
£,{t) is a centered Gaussian process such that 

\im -^E[^{t)f (4) 
t^oo log(t) ^ ' 

that is, the width of the interface fiuctuations increases in time as log(i). . 

However, a number of natural phenomena cannot be described by simple 
diffusion; e.g., the way some proteins diffuse across cell membranes or the motion 
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of a particle in systems with geometric constraints, for example on the surface 
of a perfect crystal. Therefore, it can be of interest considering systems of SDEs 
analogous to ([IJ and in which the Brownian particle interacts with anomalously 
diffusing fields. The present paper is devoted to extending the results obtained 
in [I] for system ([T|), to the case in which the interface fluctuations are due to 
interactions with anomalously diffusing fields. In other words, we will study the 
long time behavior of a Brownian particle coupled with an anomalously diffusing 
environment (see systems pT|) and (HI])). 

Anomalous diffusion processes are characterized by a mean square displace- 
ment which, instead of growing linearly in time, grows like t'^"', 7 > 0, 7 7^ ^. 
When < 7 < i the process is subdiffusive, when 7 > i it is superdiffusive. 
Diffusion phenomena can be described at the microscopic level by BM and 
macroscopically by the heat equation, i.e. the parabolic problem associated 
with the Laplacian operator; the link between the two descriptions is, roughly 
speaking, the fact that the fundamental solution to the diffusion equation is the 
probability density associated with BM. 

A similar picture can be obtained for anomalous diffusion. The main difference 
is that in nature a variety of anomalous diffusion phenomena can be observed 
and the question is how to characterize them from both the analytical and the 
statistical point of view. It has been shown that the microscopical (proba- 
bilistic) approach can be understood in the context of continuous time random 
walks (CTRW) and, in this framework, a process is uniquely determined once 
the probability density to move at distance r in time t is known ([3]- [7], [18], 
[19] and references therein). The analytical approach is based on the theory of 
fractional differentiation operators, where the derivative can be fractional either 
in time or in space (see [8] -[10], [17] and references therein). 

For /(s) regular enough (e.g. / € C{0,t] with an integrable singularity at 
s = 0), let us introduce the Riemann-Liouville fractional derivative, 

and the Riemann-Liouville fractional integral, 

where T is the Euler Gamma function ([10]). Appendix [B] contains a motivation 
for introducing such operators. For ^ < 7 < 1 let us also introduce the frac- 
tional Laplacian A^'^\ defined through its Fourier transform: if the Laplacian 
corresponds, in Fourier space, to a multiplication by — fc^, the fractional Lapla- 
cian corresponds to a multiplication by — | A; |~. ([5]) and ([S]) can be defined 
in a more general way (see [10]), but to our purposes the above definition is 
sufficient. Furthermore, notice that the operators in ([5| and ([6]) are fractional 
in time, whereas the fractional Laplacian is fractional in space. 
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Let us now consider the function p'~'{t,x), solution to 



""""•^'^rxbsi'*!^ whe„0<7<i. (7) 

W--) = f(2T^/'''»|^ -'■e„l<,<l. (8) 

It has been shown (see 21 UHl ^.nd references therein) that such a kernel is 
the asymptotic of the probability density of a CTRW run by a particle either 
moving at constant velocity between stopping points or instantaneously jumping 
between halt points, where it waits a random time before jumping again. On 
the other hand, a classic result states that the Fourier transform of the solution 
p'^{t^ x) to 

9tp^(i,a;) = iA(V(t,a;), ^ < 7 < 1, (9) 

is, for 7 > ^, the characteristic function of a (stable) process whose first moment 
is divergent when 7 > 1 (see [13); this justifies the choice i < 7 < 1 in equation 
Processes of this kind are particular CTRWs, the well known Levy flights; 
in this case large jumps are allowed with non negligible probability and this 
results in the process having divergent second moment. 

We will use the notation p'^{t,x) = p2{x) to indicate the solution to either ([7]), 
([H or (|9]), as in the proofs we use only the properties that these kernels have in 
common. 

The above described framework is analogous to the one of Einstein diffusion: for 
subdiffusion and Riemann-type superdifFusion the statistical description is given 
by CTRWs, whose (asymptotical) density is the fundamental solution of the 
evolution equation associated with the operators of fractional differentiation and 
integration, i.e. ([7]) and ([5]), respectively (see Appendix B). For the Levy- type 
superdiffusion, the statistical point of view is given by Levy flights, whose pro- 
bability density evolves in time according to the evolution equation associated 
with the fractional Laplacian, i.e. © (see [T7j). 

In view of the above considerations, we introduce the following three systems 
of Ito-SDEs: 



Id/"* d'^h'^is x) ^^^^ 



dX-^it) = X'^dwit) + Xy-\ipx-,(t),h'^it))dt 



dh''{t) = 



1 



r(27-i) 



ds 



dlh-^js.x) 
[t - s)2-27 



(11) 



dt - (px-,{t)dX^(t), 
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and 



(12) 

[ dW{t) ^ \/^^''^h^{t)dt~Lpxi(t)dX-<{t). 

Roughly speaking, the first two systems are obtained from ([l} , by replacing the 
Laplacian of the field h{t, x) in equation ([IJ2, with the fractional derivative and 
fractional integral of A/i(t, x), respectively (see (O and (|H])). The last system is 
obtained by replacing the Laplacian with the fractional Laplacian (see ©)■ In 
this way we model our anomalously diffusing fields. 

Again, 'w{t) is a one dimensional BM, <p(x) is a function in the Schwartz class 
and (px'i{t) = fi^ ~ -^'^{'t))- A more detailed motivation for introducing the 
above systems can be found in Appendix B. 

We shall denote by (t) the solution to either of the three above systems (the 
reason for adopting this notation, which might at first seem confusing, will be 
apparent in few lines). For A e (0, 1), let us introduce the scaled variables 

h^^-''\t,x) := ^h-^ (x\-^,t\-^^ , 

(^(^)(x) :=i(^(xA-i). 

For the function only, we use the convention ipa{x) := ^p{x — a), a S R and we 
set 

ip'^^\x) = \^ (x\-^ - X(^'^)(t)) ; (13) 

the notation for (pj^-* should include the superscript which we omit. 
Let also £,'^{t) be the solution to the integral equation 

c{t) = h{t) - f dsp]_MCis), r(o) = 0, < 7 < 1, (14) 

Jo 

where b{t) = w(tX^^). Notice that, in virtue of the scaling property of 
Brownian motion, the dependence of ^'^{t) on A through b(t) is only formal. 
The main result presented in this paper is a scaling limit (in fact, three scaling 
limits) of X'^-'^-'^^t) to S,''it). Also, the solution to ([Ti)) is unique by basic facts on 
the theory of Volterra integral equations, which we shall recall at the beginning 
of Section O 

Theorem 1 (first version). With the notation introduced above, we have that 
V7 e (0, 1) and ViV £ [1, 00) there exists t ~ t{N, 7) > such that 

limE sup iX'-^-^^t)- = 0, 

t<T\ logAI^^ 

where C(7) is a positive constant, with C(l/2) = 1. 
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The fact that C(l/2) = 1 is coherent with ([3]). In Section!?] we prove an 
equivalent version of Theorem [Tl namely Theorem [3l where the constant C(7) 
is made explicit. Theorem [T] says that the asymptotic behavior of X^^''*\t), the 
rescaled solution to either one of the systems ((TU)) . pT|) and (IT^ . is described 
by the function £,'^{t). Hence, we need to determine the behavior of £,''{t) for 
large t, and this is the content of the following Theorem [5J 

Theorem 2. For 7 € (0, £,''{t) is a centered Gaussian process s.t. 



lim E[Cit)f ^ const. (15) 



Forj^ (i'l); we prove an invariance principle for S^"^ (t) . LetQ{t) — e^^^ S^"' {e^^t); 
then, as e — > 0, ^2 converges weakly in C(]R+) to a mean zero Gaussian process, 
Z{t), whose covariance function is 

E{Z{s)Z{t)) = ^ y;! / du- 



T^'^c{-if Jo {t-uy-'is 



U-7 ■ 



Intuitively, this means that in the case in which the particle interacts with 
a subdiffusive field, the feedback force exerted by the field keeps the process 
localized. On the other hand, the superdifFusive field (no matter which one of 
the two we consider) is not strong enough to contrast the effect of the Brownian 
nature of the particle and the width of the fluctuation increases in time as t'^^~^ . 
Notice also that the CTRW associated with the operators and are non- 
Markovian whereas Levy processes are Markovian processes; nevertheless the 
limiting process ([T4)) is non-Markovian for any 7 G (0, 1): in the case of Levy- 
type superdiffusion there is loss of Markovianity. 

The paper is organized as follows: in Section [21 after establishing the notation, 
we state a second (equivalent) version of Theorem [TJ This version is the one 
that we shall actually prove in Section S) Section [3] contains all the technical 
estimates used in Section S) This proof is a generalization of the one used in 
[1] in order to prove Q . Section [5] is devoted to the proof of Theorem [2l which 
relies on the use of Tauberian Theorems. Finally, Appendix A provides a sketch 
of the proof of existence, uniqueness and continuity of the paths of the solution 
to (fTO|l . pT|) and (fT2|l . Appendix B contains a more detailed motivation for the 
introduction of the operators of fractional differentiation and integration. 



2 Notation and Results 

The kernels in ([7]) and ([8]) can be explicitly written both in integral form (see 
Appendix [B|) 

I pc+ico g-|a:|z^ 

p'^it, x) ^ — dz e^* ^_ Vc> and < 7 < 1 (16) 
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and as a series, 

M{z,j) 



1 

oo 

E 

k=0 



{-if 



< 7 < 1, where 



/c!r(-7(fc + l) + l)' 



(17) 



The asymptotic behavior of the Mainardi function M{z, 7) as z — > +00 is known, 

M(z, 7) ~ A(7) e-^(^) , 

with A and B constants depending on 7 ([9]); hence x) has finite moments 
of any orders, which is 

dxp^{t,x)\x\''<oo, WneN. 

We remark that this property holds when p^(t,x) is the fundamental solution 
of either ([7]) or ([5]). On the other hand, the fundamental solution of Q, namely 



Ulfcl ~' ikx 



e'^'-'dk, 7 e (1/2,1), 



(18) 



has finite first moment but divergent second moment. 

We want to remark that in order to prove Theorem [T] (i.e. Theorem [3]), we 
basically use only the following elementary properties enjoyed by both (jl6p and 



• scaling property: 
from which, setting 



p^(^,x)4p^(i,^). 



c{l) :=p7(0), 

pJiO) c(7) 



{t-sp (i-s)T' 
• there exists a generic constant C — C{'y) > such that 



c(7) 



< 



C\zfB V/3e(0,l], 



and 



c(7) 



- 1 



< 



cE 



(19) 

(20) 
(21) 

(22) 
(23) 



^This inequality can be deduced by using 11711 when referring to Riemann-type anomalous 
diffusion, see also footnote [S] When pj^ is the kernel in 1181 1. see footnote |3] 

■^The constant that appears in this inequality is equal to 1 when is either the Levy- type 
kernel or the subdiffusive kernel and it depends on 7 otherwise; see again footnote |3] 



7 



For f,gGL'^{[0,t\),f*g denotes the Volterra convolution, namely 
(/*5)(<):= f dsf{t-s)g{s). 



For TO G N, m > 2, /*(™) = /*/('""^) is the convohition of / with itself (m- 1) 
times, where we define f*^^\t) := f{t). Set K^{t) := p]{Q) and notice that 

(24) 

If we iterate n times, we end up with 

K;("+i)(i - s) := ds'p2_Amt^{s' - s) 

= /c(„)(7)(i-sr-("+^)\ n>l, (25) 

where 

^,..|(7):=chr V((n+'l)(C)) - <^^' 

To obtain the prcvioTis equality we used the fact that the Beta function B{z, w) 
can be expressed in terms of the Euler Gamma function in the following way: 

r(^)r(u;) 

T{z + w) 



Jo 



B{z,w) "= / dss'-\l - s)"'-^ = ^^.^^ """^"".^ Re{z), Re{w) > 0. 



In the same way, by setting 

P(^'^):=(^P,P7_,^W) = P*«, (27) 
(on the RHS we drop the superscript ^^'^^ for notational convenience) we have 

.,7) 



and for n > 1 



We further introduce 



J s 

P<r'^:= f ds'P^^r^P:^:^. (28) 

J s 

:= f db{s')p]_,,^^J^^), 

J s 

F^'''^\t) := - f dsKi^\ (29) 
Jo 

f -)(i) := r 

Jo 



Fr "it) := / dsPt^^A^Ki^^^\ 
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Jo 

* (A,7)p*(2) 

S 1 



Jo 



and in general 



Fi^^^\t) :^ (-1)"+! / d.xi'o^)p,:i"', n > 1 (30) 
Jo 

= - /*rf.Pi,^-'')Ffj7)(,s) n>2. (31) 

Jo 

Via the Duhamel principle (see Lemma [5]), systems ((TU)) . pip and ([T^ can be 
expressed in integral form by a unique system, that is: 

X(^'^)(t) = 6(t)+ [' ds{v(^\h^^'^\s)) 



(32) 



where 7 e (0, 1); in the above system p1{x) = p^{t,x) is either for 7 e (0, 1) 
or for 7 e (1/2, 1). For any / in the Schwartz class, (p7/) (^) ^ convolu- 
tion in the space variable. Namely, pj^^ipi'^ = (^pj^^ipi'^^ [x) — J^dyp]_g{x — 

y)tpi^\y). The initial conditions for ^ are X(^''^)(0) = /i(^'T)(0) = 0. In 
Appendix A we prove that ([32l) admits a unique solution in C(M+;M x L^(M)). 
Notice as well that from (l32t one has 



h('''r\t) = - / dX('--'\s)pU^i'\ (33) 



Following [1] page 10, we formally iterate once both the equation for X'-^''^^ and 
the one for , ([32l) i and (|T4|) . respectively. Plugging ([32l) 9 into (|32|) i and using 
we get 

x^'^'Ht) ^ bit)- fdsKi^-'h [ds{^i'\ [ds'{^i^\ [dxl':''\s")pl_^„^i^}}p:_^,^ 



'(1) y^^l ~ "y'')^J^''^^s,0 '^J^^^^ Jo ' Jo )r-s'-s"-rs" li-s-s' 

= bit) + F^''''\t) +[dsrds'p^'f\^^^\ fdx[^:^\s")pl_,J^}), 

Jo Jo Jo 

(34) 

where the subscript (1) is to recall that we are considering the first iteration of 
dSlDi. Setting Yl^^^\t) = X['^:^\t) - bit) - F^^'''\t), Y^-^^'^^t) solves 

r((f)(t) = J^dsF^'^''\s) + J^dsFt-''\s) 

+ f ds f d^P^;'^(^^, fdYl^-'\^')fil,_„^^); (35) 
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observing that Y^^s^''^\t) is a.s. in C^(R), we can rewrite the previous expression 
fory/,^'^^(i) as 



(1) 
hence 



X 



(1) 



(37) 



On the other hand, iterating the equation for and using we get 



We can repeat the same procedure n times; for n > 2 we then have: 



Xl^^-^^t) ■.= b{t)+F^'-'^\t) + ds[Fl''^^ + --- + Fi^^'r)]is) + Yl^f\t), (39) 
where 

Yl^Y\t) ■■= i-^r^' tds{^i'\ f dx\^^\u)p:_^^^^^) fds'P<-J. (40) 



Y^^^'^\t) solves the equation 



ds 



(^) 



so by differentiating, using the definition of P/^'''^'' and ([28l) . we get 



F^^'^^ + --- + F 



(A,7) 
2n 



it) 



+ (-1)"+! / dsYi:fHs)p: 



Define Aj^.^{t) as 



(42) 



n 

^(n) W T.i~^y * ^) W n>l, 0<J< (43) 

i/=0 " 

where ^K^^^-* * (t) is only formal and we set it to be equal to b{t). Then, at 

the n— th iteration of the equation for the limiting process (."'{t), we find that 
Vn > 1, 



fin)it)=AUt) + i~^r^' I d<,(.s)K*("+l)(<-.) 



(44) 
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When we write Xj^^l^^' , we refer to the expression ([39]) if n > 2 and to p7|) if 

n = 1. As for yI'^^'^'' and y(^„^)'''\ expressions gH) and ([42]) coincide with ((35l) 

and (15^ respectively, when n ~ 1. So y^'^-,'''^'' and y^^^f^'''^'' are defined by (I4ip and 
([42]) Vn > 1. 

To prove convergence of X^'*'''''^ to ^'^ we prove convergence of the n-th it- 
erates. More precisely, we prove that Vn > 1, x'^^^'^'^ converges to ^^^^ (in the 

sense of Theorem [3] below) when 7 € ^0, • 

The reason why we consider successive iterates of the equation for X^'*'''*'' 
(and hence for ^'^) is to gain integrability and some sort of regularity. Notice 
indeed that db{s)p]_g{0) is not well defined for 7 > 1/2, whereas Vn > 1 

f d6(s)K!("+i)(i - s) is well defined for 7 e f 0, . (45) 

Jo ^ \ 2{n + l)J 

Vn > 1, we further restrict the range of 7 to 7 G ^0, -^^^ in view of ([^5)) (see 
Remark li?^ and pUS)) . as well). 

Theorem 3 (id est, second version of Theorem[T]). With the notation introduced 
above, we have that V7 G ^0, and yN G [1, 00), 3t = t{N, 7) > s.t. 

lim £; s up ^ I X[^f\t) - e^^^it) r= 0. (46) 

t<r|logA| 

11 • lip, p > 1, indicates the usual LP{R,dx) norm and {pi f){x) = f dyp]{x — 
y)f{y) is a convolution in space. Now a few observations: Wt > and Vn > 1 

^i'^ = (xA-i - x(^--'Ht)) = (xA-i - x;„y)(i)) , 7 e (0, 1); 

(47) 

so actually the notation for ip[''^\ as well as the one for Kj.^j'''' and ri'^''''''^'', the 
latter defined in ([79)) . should explicitly show the "dependence" on n, but we 
omit it. This also explains why in some estimates (for example (j78p ). n appears 
on the right hand side but not on the left hand side. 
Vp > 1 there exists a positive constant C — C{p) s.t. 

||^i'^|jp<CA^-^ (48) 

Moreover, Vi > 0, 

pHt,x)<B{-f)p'>it,0), (49) 

where B{'j) = 1 if p''' is either the subdiffusive kernel or ([TSl) . and it is a positive 
constant actually depending on 7 in the case of Riemann-superdiffusion. [f| (j49p 

more detailed account and helpful plots of the kernels can be found in j9j on page 
1473; see also 1101 . As for the kernel in 1 II8II . we recall that both and its first derivative in 
space belong to L'^(R) n L°°{R), Vt > and we refer to [Hi- 



ll 



implies that 
and 



4"'"^ < s(7W-.(o), yo<.s<t, 



From (l50l). we also have 



(50) 
(51) 

(52) 



where C > is a generic constant depending on 7. 

3 Technical Lemmata 



Throughout the following Lemma we will make extensive use of the Volterra 
convolution introduced in Section[21 Notice that this convolution is commutative 
and that it enjoys the property 



(t)- / du{f*g){u), 



duf{u) * g 



which easily follows after a change of variable. Indeed 

r dufiu)^ *g {t) = Jds (^J^ dufiu)^ g{s) 

t nt nt pV 

dsg{s) / dvf{v ~ s) ^ dv dsf{v - s)g{s) 
Js Jo Jo 

t 

dv{f *g){v). 



Lemma 1. For n d N, n > 1, consider the integral equation 
/i(„)(t) - (]K;("+i) * /^(„))(i) = g{t), g e L\[0,t]), 7 G (^0, 



(53) 



(54) 



Call >cj^^^(t) the solution to when the forcing g(t) is taken to be equal to 
"^(Tt)(^) ^ ^^([Oj^]) O'lT-d (i) the solution to the same equation with a different 
forcing, say ^^^^ (t) . Namely: 

<)(0 + (-1)" 1^ ds^l^{s)k^„^{^){t - .)"-("+i)7 = Al^^it) (55) 



and 



^(l)(i) + (-l)"y^ ds,l^is)k(^^^i^)it - s) 



ri-(n+l)7 _ 
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If the two forcings (t) and Gj^^^ (t) are related through 

(AJ,^^ * K;("+i))(<) = (-1)"+! dsGl^^is), (57) 
then the same relation holds true between the corresponding solutions, i.e.: 



The proof of this lemma is an immediate consequence of some basic facts in 
the theory of Volterra integral equations, which we recall here. For more details 
on this theory we refer the reader to p^]. For some T > 0, let g{t),IC{t) G 
L^([0,T]). Then the solution h{t) to the equation 

h{t) - / dslC{t - s)h{s) = g{s) 
Jq 

exists and is unique and can be expressed as 

h(t)^g(t)- f dsH{t-s)g{s), (59) 
Jo 

where 

oo 

H{t-s) = - 

i/=0 

When the kernel IC{t) is not in L^, the solution to ([59l) still exist and is unique 
provided that for some rt e N the iterated kernel /C*(") is bounded on [0,r]. 
The proof of this fact can be found in [TT], Section 1-12, where kernels of the 
form IC{t) = t", with a G (0, 1) are considered. 

Proof of LemmaUi For 7 g (0, n/n + 1), the kernel of equations (|55l) and ([56]) is 
a bounded continuous function, so the standard theory for kernels in applies. 
Thanks to together with dSS]), (HH) and (HZ]), proving ^ boils down to 
proving 

(-1)"+! 1^ ds (h * Gi^^) is) = (k;("+i) * H * A^^) it). 

Such an equality holds true because, by the commutativity of the Volterra con- 
volution, the right hand side is equal to 



* AJ^^)] (t) = H{t - s) (k;("+i) * Al^) {s) 

= (-l)"+i f dsH{t~s) f gl){s')ds'; 
Jo Jo 

now the conclusion follows from property ([55)) . □ 



13 



In the following lemma and throughout the paper we will be using the no- 
tation = f{k) and C{g{t)}{^) = g'^i^J.) for the Fourier and the 
Laplace transform respectively and we will superscript "for the Fourier-Laplace 
transform. 

Lemma 2. For < 7 < ^, let v^{t,x) he a solution to 

a..^(M)==^^y/.(^^^ + F(,x) (0,00) xM 

wT(0,x) = ^^(a;) {0}xR 
and, for ^ < •-f < 1, let it be a solution to 

dtv^t, x) = ds + F{t, x) (0,oo)xR 

v^{Q,x)^vl{x), {0}xR 
where w^(a;) G C(M), F{t,x) G C(R+ x K). Then 

v^{t,x)^ / dyp'^{t,x - y)vl{y) + j ds I dyp"^ {t ~ s,x ~ y)F{s,y), 



with p^it^x) the kernel defined in MOjl . 

Proof. Let us observe that Duhamel principle for the heath equation (i.e. the 
parabolic equation associated with the Laplacian) can be expressed as follows: 
if u{t, x) is a classical solution to 

dtu{t,x)^\l^^u{t,x)+F{t,x) (0,00) xR, F eC{M.+ y.^), 
u{Q,x)=uo{x) {0} X M, uoeC(M), 

then its Fourier-Laplace transform satisfies 

u{Q,k) + F{^i,k) 

u{n,k) = , (60) 

/i -j- 2^ 

where (/i-f A;^/2)^^ is the Fourier-Laplace transform of the fundamental solution 
of the diffusion equation, i.e. of the heat kernel. 

Now let us recall that the Fourier-Laplace transform of p'^{t,x) is given by 
(|134p (in p34p take ci = 1, see Appendix IB|): also, p}~'^^ v^if^, k) is the Laplace 
transform of {v"' {■ , k)) when < 7 < i, whereas for ^ < 7 < 1 it is the 
Laplace transform of If{v''{-,k)) (see Appendix [B|) . Hence 



C{dtv''{t,k))^ dte-^''dtv{t,k) 
Jo 

-v{0,k) + pv'^{p,k) = -cik^fi^~-^'^v{n,k) + F{n,k) 
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z)(0,fc)+f(M,fc) 

which is precisely what we where looking for (see (|134p and □ 

Lemma 3. ViV > 1 and < 7 < 1, let p, q and r he real numbers greater than 
1 s.t. + = 1, 9 > max{iV,r} and r^^ — q^^ < {2j)~^ . Let v{-) be an 
J^s-adapted process in C(R+, L''(IR)) and ip a random variable in L^(IR). Then 
there exists a positive constant C = C{q,r,^) such that 



N\ W 

E / db{s)p2^_As)) 



<C{t2-t,Y{EUf^y (e sup \\v{s)\A 

\ tt<S<t2 J 



for any ti < t2, where /3 = and = 5 - 7 (7 ~ ^) • 

Proof, (sketch) The proof is identical to the proof of Lemma 3.1 in [T], so we will 
not repeat it. The additional condition r~^ — q~^ < {2j)~^ is an integrability 
condition and comes from the fact that 

\s\\pl_J\l,<oo^l-U^, (61) 

where r' is such that p- + 7 = 1 + | (see page 12 in ^Ij). (161]) follows from the 
scaling property (jl9p in the following way: 



t2 



ds\\pl^J\l, = / ' ( I dypU.yKt2 .s)- 

Jti 1^2 — S) ' \J]g^ 

C I ' ds{t2 - s)2^(7^-i) < 00 ^ 27(^ - 1) > -1 

111 

r (7 27 



□ 



Remark 3.1. The extra condition 7 ~ ^ < ^ *s automatically satisfied when 
7 G (0, 1/2]. It is non empty only when 7 £ (1/2, 1). 

In the remainder of this section and in the proof of Theorem [3] we will very 
often make use of the following simple observation (sometimes without saying 
it explicitly). 

Note 3.1. Let {^,p), {ft',p') two (finite dimensional) measure spaces, f : il x 
f2' — >■ R o positive function and m a real number greater or equal to 1. Suppose 

F{y) := / diJ,{x)f{x,y) < 00 for a.e. y € fl and 
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d^J''iy){ / dfi{x)f{x,y)] < oo. 



The 



d^Ji'{y) 



F{yr-' / dti{x)J{x,y) 



dfi{x) / df/{y)F{yr-'fix,y) 
n Jn' 



< I I d^,'{y)F{yy 
In' 



dfi{x) / dfi'{y)f{x,yy 
n \Jn' 



having applied Holder's inequality with m/ {m — 1) and m. Looking at the first 

and the last line of the above equations and dividing both sides by [J^^, (J^^ /)"'] ™ 
we obtain 



dfi'iy) / dpi{x)f{x,y)] < / d^^{x) / dfi' {y)f{x,yy 
n' \Jn / \Jn \Jn' 

(62) 

When (fl,fi), (fl',fi') are just M equipped with the Lebesgue measure, the above 
inequality reads 



dy / dxf{x, y)] < { dx( dyf{x, y) 



If instead (f2',/^') is a probability space and (JlyfJ.) is the time interval [0,T] with 
the Lebesgue measure, inequality W2jl implies that Vi G [0,T] and N > 1, we 
have 







N 




E sup 


f dsf{s) 


< E 


rds\f{s)\ 


t6[0,T] 


Jo 




Jo 



N 



<T^ sup E\f{s)\^. (63) 

sG[0,T] 

In the remainder of this Section, C is a constant that does not depend on A 
or 6, although it might depend on a positive power of T. Also, in the proofs we 
assume for simplicity T > 1, even though all the results are true for any T > 0, 
hence they are stated in such generality. Even if we assumed T > 1, this would 
not be restrictive in view of the fact that the result we are concerned with is a 
long time result, more specifically T ^\ log A | with A — > 0. The case 7 = 1/2 is 
not explicitly considered in Lemma |4] and Lemma [S] 

Lemma 4. VA^ > 1, < 7 < 1 and ( e ^0, 5^^, there exists C > such that: 

< CT'^'^A^-^-'^, r>0,Ae(0,l). (64) 



sup I E 

0<s<t<T 



K. 
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Also, Vn > 1, iV, 7, ^ as above 



E sup 

V te[o,T] 



N 



T > 0, A G (0, 1). Moreover, for the displacement of the center we find: 

( sup E sup |x;„Y^(i')-x;„Y^wrV 



yte[o,T],t+T<T t'e[t,t+T] 
< C (t^ + 



r,Ae (0,1) andT>0. 



(65) 



(66) 



Proof, (sketch) follows from Lemma [3] and where in Lemma [3] we have 
chosen ^ — | = CS ^0,^^. Having in mind Note 13.11 from ([^ and 

((M)) . using we have 



i? sup 
, te[o,T] 



AT 



From (1201), dMl) and dsn we get 



(67) 



sup E 

\te[Q,T] 

so that, again by (|63 



< T^^r"-"''A^~'^-\ n>l, (68) 



i? sup 

, te[o,T] 



< c'T^^+ir"-"''A57-';-\ (69) 



Also, from (gS]) and (|52|), 



sup 



sup ( i? 

te[o,T] 



Af\ TV 



+ sup ( E 

te[o,T] 



1 



+ CT"-^("+i) / dt sup E 
Jo \se[Q,t] 



N 



By the Gronwall Lemma and we then obtain that Vn > 1 



sup E 
>ie[o,T] 



(70) 
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hence 



E sup 

, te[o,T] 



N\ jv 



(71) 



When n — 1 ([^5)1 is a straightforward consequence of ([57|) . (p7|) . ([7T|) and the 
fact that 

£ sup I 6(t) r< CT^. (72) 
te[o,T] 

When n > 1, we first rewrite ([M)) as follows 

(t) = 6(0 + ^^0^''"^ it) + ds [f^''-'^ + . . . + ^(a,7)] (,) + ^* dsYl^^f) (,) , 

(73) 

and then ((65|) follows from (|67|) . (|72|) . (|69l) and ([71]) . By acting in a similar way 
we find the following estimates: 



sup E sup 

.t<£[0,T],t+T<T t'e[t,t+T] 



N 



sup E sup 

tG[0,T],t+T<T t'e[t,t+T] 



sup E sup 

te[0,T],t+T<T t'<£[t,t+T] 



N\ N 



So, recalling that for the BM b{t) 



E sup \b{t')-b{t)\^<CT^, 
f'e[t,t+T] 



66l) follows. 



(74) 



□ 



Lemma 5. ViV, n > 1, < 7 < 1, C G (^0, , T > 0, A, (5 G (0, 1) i/iere exists 
a constant C > such that 









\E sup 


I s.s — o 






Js 





\E sup 



dsp2_MK, 



s.s — 6 



<c 



(75) 

^^(7T+l)(l-7) 

(76) 
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E sup 

V t6[0,T] 



t N\ — 



< 



CX^-^C-\CTi-^'>''-^\ ^^^^ 



Moreover, VAf > 0, we have 



sup [E 

se[s,T] 



<C 



CT 



where 1 is the indicator function and 

s-S 



i-5 



e 
(78) 

(79) 



Proof. The proof of the bounds (|75j), (|76l), (|27l) and (Hlj) is done by foUowing 
P, pages 16-18, so it shall not be very detailed. Recalling (|47)) . we have that 
Vn > 1, 7 e (0, 1) and for < s < Kf;'^' can be expressed as 



(„ + 1)(1-t) 



(A,7) _ 



if' 



^(^), / db{.s')pt,,^'^'>) 



+ {v['\j^db{s')p1_,{^\V~v['^)). 
Observe also that Va, 6 G K and Vm > 1 

ll-^b - VaWm < ll'^'llm \ b- a \ , Lph-^ " /l)- 



(80) 
(81) 
(82) 

(83) 



Let us start with proving (|77p . We decompose if^'^'''''' according to the prescrip- 
tion (|80p - ([5^ : recalling the notation (IT5|) . the term coming from (|5T|) becomes 



Using Lemma 131 we have 



E 



(^w-^(^), / db{s')p:_,v^^)) 



N\ — 



(A) I 



with r,p and v to be chosen according to Lemma [H By (|83|) . (|48|) . and (|65l) . we 
obtain that 



sup \E 

sG[0,t]\ 
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having chosen + = ^ — C- For v' ^^nd 5' such that + l/g' = 1, we 



have 



r p 
t 







< C 



ds 



- s)^''' 



<3 \ 9' 



SO that 



sup 

, te[o,T] 



Jo Jo 











< c 


if 










Jo 











1' 




/ ds 


(^(^)-^(^), r db{s')p^^^,^(^^) 






Jo 


Jo 





N\ Ng 



< C sup \E 



(^i^)-^(^), / db{s')ps-s'V^'^) 



Nq'\ Nq' 



(„ + l)(l-T.) 



The addends ([80)1 and ([82|) can be examined in the same way, so we leave it to 
the reader. We now very briefly show how to obtain ([75)) . We decompose again 
according to (l5ni) - ([5^ . For the term coming from ([5(I| . by exchanging 
the order of integration (which is now allowed) and integrating by parts, we get 

15 Js-S 

For the term coming from (j8ip . we have 

N\ w 



E sup 

, *6[<5,T] 



E sup 

, te[5,T] 



-5 



having applied Lemma [3] with the choice 7 — ^ = ^ — d C ^ (o, ^ j , and 
as well. In an analogous way, for the term coming from ([82]) we obtain 



E sup 
, tels.T] 



ds( 



N\ — 
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(|76p results from applying again the same technique so we won't present the 
proof. 

In order to prove ([78|. let us express ri^''^'''^'' as 



where 



r(A,<5,7) 



{x) 



t-5 



) / dyipiy) {pt^ [X{x - y - x[^:^''\s))\ - p]_m} • 



By a change of variables and using the scaling property (jl9p , we can rewrite 



t-S 



where 0(7) is defined in ([2D]). We now use the bounds (|22p and ([25)1 . More 
precisely, setting z = A — — -'^(^j''^ (s)^ /{t — s)'^, we estimate the integrand 
above in the following way : 



c(7) 



p7(£) _ I 
c{j) 



-I <C when | x \> X-^/^ 
<C\z\ when I X \< A"^/®. 



So, following jlj, pages 15-16, we apply the Burkholder inequality (tl2j) and we 
get 



i?i/t(x)r<ci 



{\x\>\-^/<^} 



ds 



(t's)'' 



+ Cl{|:r|<A-l/8}^ 



t-S 



ds 







- ^l{|2:|>A-i/s} |l{o<-,<i/2}*^ ^'^ + 1{1/2<^<1} 



^1-27 1 



-'V s6[0,T] 



(") 



t-5 



(t-s)'-^ 



- ^l{|2:|>A-i/s} |l{o<^<i/2}*"^ ^'^ + 1{1/2<^<1}'^^ ^^1 
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+^l{,.,<.-vs}A- (a--/« + 1 + A^-^-) e^-'"— 

where in the last inequality we have used (|55|) . If we choose C = 1/8 in the 
above, we obtain 



E 



{|x|<A-l/S} 



N\ -N 



<C\M^\El dx 
Moreover, for any M > wc have 



<A-l/8} 



(A,<5,7) 



N 



,CT 



(„ + l)(l-T.) 







E 




V 





N\ — 



< x~ 



2N-1 



E 



>A-l/8} 



E 



dx 



r(A,5,7) 



{x) 



2N\ 2N 



y ■'{|x|>A-l/8} 



1+X^ 



- l-^{0<T<l/2}'' -^{1/2<^<1}" j 



This concludes the proof of ([78)) . □ 
Lemma 6. VO < s < A e (0, 1), /3 e (0, 1], n, TV > 1 and e (0, 1) we /law 



e( sup (i-s)(i+«^|p(^'^'-p7-.(0)| 

Vo<s<t<T 



^feo, V(5 G (0, 1) anrf for any Q > 0, we have 



(84) 



^tSTr] ilo Pl-sm] I < C (X^-'^ + XQ) 



(85) 
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Sketch of proof. Using the definition of p/^''''-' (|27p . by change of variables and 
the scahng property ([T9| . we have 



<p7-s(0) JJdxdy (p{x)ip{y) 
From (HI), then 



c(7) 



A(x-y+x;:;'w-<7\^))' 



(n) 



(86) 



We now want to use ([M)) in order to conclude; though, ([M]) holds only for > 1 
whereas /3 is in the range fi € (0, 1]. We don't want to choose fi — 1 (see (|109p 
and comments after it), hence we first apply Young inequality with p — 1/ 13 and 
get 



(N + |y| + KV^W-^£r^(s)| + i) X' 



{t - s)^P 

and now (|84|) is a straightforward consequence of (|66p . To get (|85|) . we use 
again the bounds ([^^ and (P^. this time in the following way: setting z = 

X(x-y + V) - lit - we estimate 



pM. _ 1 

c(7) ^ 



1 < C when | a: |> A'^ 
<C \z\ when I x |< A-^ 



So, from 



we have 



t-s 



< C 



t-5 



ds 



[t - S)27 



t-<5 



- S)7 



< CI, 



t-5 



ds 



{|x|<A-i} 7p (t-s)27 



+ C1, 



ds 



{Ni>A-i} (t-s)T 

now follows from (p5|) . 



(c+|x;„y)(o-4^-^(s)|) 



tpix)\x 



□ 
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4 Proof of Theorem [3] 

We recall that C is a positive constant that does not depend on A and S, though 
it might depend on a positive power of T. Also, for simplicity, all the proofs are 
presented for T > 1, even though the statements are clearly still valid for any 
T > 0. Since it has already been treated in [T], the case 7 = 1/2 is not explicitly 
considered. 

The intuitive idea that motivates the structure of the proof is based on the 
observation that, "morally", things go as if P/^ ''''' were converging to p7-s(0) as 
A "?> (see Lemma|6]); formally, this can be obtained by thinking that, as A — > 0, 
(p['^^ —7' So- While such an idea is not hard to turn into a rigorous argument, 
one of the main technical difficulties is encountered when trying to do the same 
thing to get an intuition on what -ft'^'o^'' ought to converge to. If in the definition 
of Xj^'''^ we replace ip^^"' with Sq and exchange the order of integration, we find 

that K^.^q''^ should converge to /J db{s)p'J_^{0). The problem is that we are not 
allowed to exchange the order of integration (see comment after (3.5) in [Tj) 
and that /J db{s)p]_g{0) is not weU defined as a process in C(K+) when 1 > \- 

In the same way, Vn > 1, -F'^'^'^'' is well defined for any 76 (0, 1), whereas the 
object it converges to is not (see (|116|) and (|45|) ). 

The proof goes as follows. Vn > 1 we introduce the process rij^-^{t), solution 
to the equation 

<) W = GUt) + (-1)"+' / dsvJ„^{sW("+'\t -s), < 7 < ^ (87) 
where 

2n „t 

Gj^^{t):^Y.'^~ir^' db{u)K;('^+'Ht-u), n>l,0<7<^. (88) 

We now observe that Lemma [T] can be applied to Cln)^ defined in pij) . and ijjn)- 
In this case the forcing terms are ^J^) and G^^^ , respectively, and we can easily 
prove that they are related through ([57| . We can in fact show that the i — th 
addend of AJ^-^ is related to the i — th addend of Gj^^^ through ([57]) : all we need 
to show is that \/v G 0, n, 

i^iy (k^"^ * & * {t) = (-i)"+i^ds(-i)''+"+i^ d6(u)K;(''+i)(s-u), 

which is a straightforward consequence of the definition of K^^™^ given in (|25p . 
together with the following equality 

(ik;("+iU&) (t) =^ ds^ d6(M)K;(")(s-M), n>l. (89) 
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Hence, Lemma [T] gives 

(-1)"+! (e;:,)*K;("+i)) {t)= fdsvi^is). 



(90) 



Recall that the definition of xj^"^^ is given by ([5^ for n > 2 and by (1571) when 



1. Using (IMl) . we look at the difference between x'^^'^'^^ and ^J^y 



3 = 1 ^•'^ 



V / dsFj^-^\s) - (-1)^+1 (K*(J+i) * 6) (t) 

.7=1 ^ ■ 



(91a) 
(91b) 

(91c) 
(91d) 
(91e) 

(91f) 



where for n = 1 the sum in (j91b[) (and in (|91ep ) is understood to be equal to 
zero. As we have already said, we want to prove that Vrt > 1, X^^^^j'''-* converges 

to ^^^j for 7 G ^0, • To this end, let us further expand the integrand in 

(|91f|) . using the fact that Y^:^^''^ solves equation (|42p : 

<^ 

(92) 



+ („i)(n+i) r dsYl^f\s) [P4*i"+'' - K*("+i)(t - s)] , (93) 



and G^^^ (t) is defined in 
Let S € (0, 1). From now on we assume t > S. 

Remark 4.1. We omit to study the case t < 5 because it can be treated in 
the same way as it is dealt with in 'll, where it is presented explicitly, see in 
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particular (3.23), (3-44) ^.i^d (3.45) in |7]/. In other words, what we actually 
show is that the estimates in MO&i) . \113\) - (TJ7^ and \121\) are valid when the 
supremum is taken over the interval [S,T] (more precisely, in the case of hll3\) - 
ill7\ ) and il21\) the supremum should be over [A", T], because at that point S will 
have been chosen to he equal to A°, see lines before jllT^ )■ Though, by acting 
as in we can show that the same estimate holds true when the supremum is 
taken over the whole interval [0, T] . Hence from now on we will assume t > S 
in order to streamline the notation and the presentation of the proof. 



We use the decomposition 



dsptsmi^^^ tdspts{0)K, 
Jo 



(A,7) 
s.O 



f-<5 



+ / dsp]_M{d'\ db{s')p:^M') + dsptMivi'K db{s')p:^M'), 

JS Jo JS Js-S 

which follows from the definition of A'^^q''', to rewrite the difference between 
f[^-''^ and /g db{s)K*^^^\t - s). For 7 e (0, 1/2), 

F[^^^\t)^ fclb{s)K<-'\t-s) 
Jo 

^d.p(^'-)<o^)- fdbis)Kf\t-s) 



N 



(94) 



<C 

+c 
+c 
+c 
+c 
<c 
+c 
+c 



t , N 



ds P, 







Pl-siO)) Ki)f^^ 



(95) 



^(A,7) 
^s,0 



f dsp1_MK, 
Jo 

dsp]_M{v>i'\ I dbis')p:_^,^i^^)^ I dsptM I db{s')pl_AQ) 



N 



s-5 



dsp1_M i db{s')p:^AO)- I dbis)K*('Ht-s) 



N 



(96) 



dsptM{^^'\ / dbis')p:_^,^i^^) 



N 



dslR 



N 



fdsp^.MKi'if^ 
Jo 



dsp]_MKi'r-\ 



N 



N 



c 



c 



dsptM^i'-'''^ 



N 



(1) 



AT 



(97) 
(98) 
(99) 
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where in the last inequahty we used the definition of F 
we set '^'l^^7\t) to be the difference in (IMl) . namely 



given in ([79)) and 



¥^l-;\t) := j\sp1^Mf^ dh{s')pl_AQ)- j\b{s)K<^Ht-sl 7 e (0,1/2). 
For n>l, we define 

*g':l)W:=^*d^p7-.(0)^S7^W, 7e(o,^). (100) 
In the same way, by using ([3T|) . ([50)) and (jlOOp . we have 

Ff'^^(t)+ f dh{s)K<'^\t - s) 



< C 



p. 



AT 



(101) 
(102) 



N 




+ C 


fdsptMf 




Jo JS 




N 



N 



dsptM I ds'p:_A0)K^'.7ls 



c 



N 



(103) 
(104) 



Remark 4.2. We will show that the terms in i98\) and the first addend in i99\) 
(hence also the addends in il03\) and the first addend in ^104^ ) are small for 
7 € (0,1) (see I1106\) . IjlOT^ and |y6] )j. The reason why we need to iterate the 
equation for X^^'"'^ and an infinite number of times comes from ^[^'^'^^ (see 
114 5\ l and il05\) ). We will in fact prove that 



E sup 

, tG[5,T] 



*&^(0 



N 



(105) 



Also, we will show that |ff7| ) is small when 7 G (0, 1/2) and ilOS^) is small for 
7 e (1/2, 1), see nJd\) and I77g)) . 

Let us now address the points mentioned in Remark 14.21 in the same order 
in which we listed them. 

For p,q > I s.t. + q^^ = I and P7 < 1, we have 



['dsp]^MK^ 
Jo 



(A,7) 
s,0 



N 



< c 



ds 



(t-sy 



Since t > 6, 



ds 



P7 



< 



(S-s) 



PI 



ds I K 



(A,7) 
s,0 
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hence 



E sup 
te[5,T] 



N 



ds I K 



(A,7) 
s,0 



< C5— 



'-N 



E 



ds I Ki)r^ r 



N\ q 



<C5^''5^ sup (e \ K^^^^ 

se[o,T] ^ ' ^ 



where in the last inequahty we used Note 13. ll If we choose p = ^^j-!- and q = 



by using ([M)) we get 



E sup 

te[o,T] 



fdsp1_MKi'^^ 
Jo 



N 



<C6^-'''X^-'^-\ 7e(0,l). (106) 



By the same sort of trick used to get (|106p . we also get 



d.s'p':_Ao)Ti^-''''^ 



Therefore, using ([75]) . we have 



N 



< C sup [E 

s'e[s,s] 



p(A,<5,7) 



Nq\ q 



E sup 

sG[A',T] 



N 



< c 



A 4 + A 2t 4 



(n + lXl-T.) 



^ -^{0<T<l/2} " -^{1/2<7<1}^ 



e--"^^"^"", 7^(0,1). (107) 



Notice that on the right hand side of the above equation, n appears because 
is contained in the definition of ri'''''^'''''*, see ((78)) . (l47t and the comment 
after it. 

As for the first term in (|99p (respectively, the first term in (|104p ). we just use 
(|76|) in Lemma [5j In order to prove (llOSp , we show in some detail how the 
estimate for '^'^^i^'' is obtained; the way one gets (|105p for n > 1 should then be 
obvious from the definition (llOOp and using (^5]) . Recalling that we are assuming 
t > 5, using (I24p and exchanging the order of integration in the definition of 

"^^^C^ we have 



(1) 



^ii)(tr-''^^ I db{s) [ ds'p]_AO)pl_M- I dh{s) [ds'p]_AO)p:,_M 
Jo Js+s Jo Js 

= f db{s) rds'ptAO)pl^siO)~ fdb{s) f ds'p1_A^)pl_M 



t-s 



(108) 



Now we can estimate the two terms in (|108p separately. In both cases we first 
make a further change of variables and then integrate by parts the stochastic 
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integral. We show how to handle the first, for the second the procedure is the 
same: 



rt-S 



s+5 













f db{s) [ 






Jo Jo 



< 



< 



b{t~5) / dup}_^{Q)pl{0) 



h{t-5) / dup}^^{Q)pl{Q) 
Jo 

Notice now that from (l24l). 



~ sup I b{s) 

se[o,t-S] 



- sup I b{s) 
se[o,t-5] 



t-5 



OS Jo 



dup]^jO)pZ{0) - / dup]_^iO)pZiO) 



dupj_jo)pZ{o) ^ cs'-'-' 



and, since t > S, 



du 



it-u)''u-' 



< C 



du 



id-u)^u'' Jo 



dupj^j0)pl{0). 



So, after dealing with the second term in (|108l) in an analogous way, (|105l) 
follows by using (17^ . 

Let us now turn to ^ and ([TU^ . Let /3 > 0; then for applying the 

Holder inequality, we have 



ds Pi 



AT 



< sup 

0<s<t 



N 



ds 



K 



(A,7) 



ds 



(109) 

Looking at the last integral in (|109p . we need to impose the integrability con- 
dition /3 < — 1 + 1/7. Taking the supremum for t e [0,T], the expectation of 
both sides, using (|64|) and (|84p . we then obtain that for 7 e (0, 1/2) and for any 
> 1, 



E sup 

te[o,T] 



ds 



p1_M\ kYo 



AT 



(„ + 1)(1-t) 



, (110) 



where we have chosen /3 — 1/2 in ([84l) . We can make such a choice for /3 because 
when we study the difference in (|94|) . and hence (|97p . we take 7 G (0, 1/2), see 
Remark 14.21 When we consider (|102p . we can't mimic what we have done for 
(^7)) : in fact from (I109P we get that the left hand side of piOp is bounded 



by A 



exp(CT("+i)(i^^)). When we impose the integrability condition 
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/3< -1 + 1/7 and/3 + i-C-2>0, ;3g (0,1], we find that these two conditions 
together cannot be satisfied for all 7 S (0, 1) (actually they hold at most for 
7 € (0, 2/3)). So, when 7 € (1/2, 1) we need to do something else. 



ds / ds'pI_..(0)(pi,t'^^-p7-.(0)) K^:o 
Jo ^ ^ 



N 



< c 



-c 



t-5 



ds ds'p:_Ao) 



p 



K 



s',0 



AT 



ds / dsVI_,,(0) 
t-s Jo 



P. 



(^,7) „i 



K 



(A,7) 



JV 



< C sup 

s6[0,T] 



ds'p:_Am^^g^ 



N 









sup 


I ds 


p.^^^-'^-ptM 


\ tG[5,T] 


Jo 





t-6 



(111) 



(112) 



where in the last inequality we have used (jSG]) and then ([23]) . By ([771) and (185 
we then have 



N\ TT 



A +0 



£; sup / ds / ds'p _ , (0) p^^'-f)- pI_ (0) ^-^'-^J \<cy 

\ te[o,T] Jo io " ° ' ' = ,0 y 

(113) 

If in (fTU5)) . (fTUSl) . ([TU71) and ^ we choose (5 = A and M > 0, recahing (ITTU)) 
we have that for 7 G (0, 1/2) and VA^ > 1, 3^(7) > s.t. 



E sup 




\ i6[0,T] 


Jo 



N\ JV 



< 



CxH-f)eCT^''-^\ (114) 



Via dSI]) and dSO]), this implies that for n > 1, 7 G (0,1/2) and ViV > 1, 
35(7) > s.t. 



E sup 

, te[o,T] 



Fi^'^Ht) - (-1)^"+'^ / d6(s)IK;("+i)(t - s) 



7V\ N 

<CA''(^)e'^^'''"^'. 

(115) 

On the other hand, if in (fT05)) . ([T06l) . (IT07| and ((761) we chose (5 = A°, with a = 
2!^, and M > ^i^^, recalling we find that Vn > 2, 1/2 < 7 < 

and iV > 1, 3/(7) > s.t. 



E sup 

I te[o,T] 



F^^^^Ht)-{-lY''+^n d6(s)IK*("+i)(t ~s) 



(116) 
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Note 4.1. We want to stress that the above estimate 11116]) is needed only for 
n >2 and 1/2 < 7 < whereas HI 15]) is valid for any n>l and 7 G (0,1/2). 



In other words we will not need an estimate on 
for ^> 1/2. 
Set now 



(A,7) 



{t)-gdh{s)Kf\t-s) 



*(0)W<''^^ / dsh{s)pJ^M- / ds / d6(s')pl.,,(0), 



then 



(A,7) 



dsh{s)pJ_M 



N 







N 


j'dsTi^^''^^ 


N 




N 




< c 


fdsKi:^ 


+ C 


+ C 




+ c 






Jo 








Js 





N 



It is easy to prove that 



E sup 

, te[5,T] 



N 



So by dMl), dZHl) and ([75]), by choosing again J = A'', a 



27-1 



{0<7<l/2} ' 27(1-7) {1/2<T<1} 



and M > • 1 



4(27-1)^ 



E sup 




V te[o,T] 


Jo 



^,,ox I -7T — ^1/,,,^ ,we get that Vn > 1, < 7 < 

{0<'y<l/2} ^^_-yJ {l/2<-y<l} ' G — ' / n+1 

N\ W 

(117) 



and ViV > 1, 3 771.(7) > s.t. 

dsbis)p]^M 
We wiU also need the foUowing estimate 



< CA' 



>(7)gCT'"+i«^-^' 



E sup 

, t6[<5,T] 



("5,7; 
(") 



(^) 



N\ N 



(118) 



This inequaUty can be worked out with calculations analogous to those needed 
to obtain (|105p . hence we omit them; roughly speaking, looking at (|105p . (jllSp 
is correct thanks to the further integration. Also, it is what one would expect in 
view of the fact that dsb{s)K*'-'^^^'> {t — s) is defined for any 7 S (0, 1), as op- 
posed to c?6(s)K*''"+^'(i — s). With this remark in mind, it is easily seen that 
with the same steps that lead to an estimate on 

using this time ()113p and ()118p . we have that Vn > 1, 7 G ^0, i;;^^ and 
3t = T{-f,N) > s.t. 

Fi^^'^Hs) - (-l)("+i) / ds6(s)K*("+i)(t - s) 



Fi^'^\t) ~ (-l)(»+i) d6(s)]K*("+i)(i - s) 
ViV > 1, 



lim E 



sup 



= 0. 

(119) 
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The last ingredient that we will need in order to conclude is the following esti- 
mate: Vn > 1, 7 G (0, and ViV > 1, 3 ^(7) > s.t. 



E sup 

, te[o,T] 



N\ TV 

(n+1) _ ^.(n+l) _ < ^;^rf(7)gCT<"+i)(--) ^ 

(120) 

which is obtained by combining ([84]) and (j71l) when n — 1; when n > 2, we act 
like in ((TTT]) - (fTT2l) and then use ([85]) and ([71]). 

From the definition of given in using pT5| . ([TTBl) and ([ng| . it is 

straightforward to see that 3c?(7) > s.t. 



E sup 

, te[o,T] 



TV 



(121) 



for any n > 1, 7 e ^0, -^^^jp^j and iV > 1. Hence, the Gronwall Lemma applied 
to (IMl), gives that Vti > 1, 7 e (o, and > 1, 3t = t{-i,N) > s.t. 



Hm E 



sup 



TV 



= 0. 



(122) 



t<r|lnA| 



Finally, looking at (|9Td| . (|9Tel) . jail, thanks to ([TT7|) . (fTTOl) and (fm]) . Theorem 
[3] is proven. 



5 Proof of Theorem [2] 

In the diffusive case, the integral equation 1^ is explicitly solvable. To our 
knowledge, ([14]) cannot be solved for 79^^. However, considering the associated 
Green function, that is, the solution of 

F^t) = 1- f dsp]_MF^is), < 7 < 1, (123) 
^0 



one gets 



db{s)F'f{t- s), 0<7<1. (124) 







Notice that the theory of Volterra integral equations for kernels with bounded 
iterates implies that the solution to (|123p is unique, as commented at the be- 
ginning of Section [31 after the statement of Lemma [T] 

Lemma 7. For any < 7 < 1 , the following holds: 

1 ™ „ . ^ sin(7r7) 
lim t^-^F^(t) = V4^, 125 

where 0(7) is defined in 
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Remark 5.1. Since c(l/2) = (27r)-i/2, Lemma\]\ is an extension of Theorem 
2.2 in IJJ. When 7 = 1/2, it provides an alternative proof of such a theorem. 

Proof of Lemma By taking the Laplace transform of (|123p we obtain that 
the Green function F'^ has Laplace transform 

(F"')*(u) = f^—r-. 7- (126) 

Provided that F^{t) is monotone decreasing, the Tauberian Theorem for densi- 
ties (see e.g. [M]) gives 

hm t'~''F''{t) = -1- ]im n^F'^)*{^i). 

Therefore the only thing we need to show is that F^{t) is monotone decreas- 
ing. We recall that a function is completely monotone if and only if its even 
derivatives are positive and the odd ones are negative. Furthermore, a function 
is the Laplace transform of a positive measure if and only if it is completely 
monotone (see again [Mj). We think of dF^it) as a (a priori signed) measure 
on M+ and introduce 

/•oo 

:= e-f^dF'^it) = 1 - ^{F^)*{pl). 

Jo 

By we have 

/.i-7 + c(7)r(l-7)' 

The function (0, 00) 3 fi — > ^^"'^ is positive and has completely monotone 
derivatives. For ^ > the function (0, 00) 3 x — > A{A -I- x)~^ is completely 
monotone. Hence (see [14]), the function $*(/i) is completely monotone and 
we are done. □ 

Proof of Theorem O By (|124p we get 





so (|15l) is straightforward. In order to prove the invariance principle in Theorem 
1, we first need to prove tightness of the process £,]{t). From (|124p and (|125p few 
computations show that for each 76 (^j l) there exists a constant C = C(7) 
such that 

c— >0 

Since ^7 is a Gaussian process, we can first obtain a bound on the higher 
moments, thus getting tightness from the Kolmogorov's criterion. Finally, the 
convergence of the finite dimensional distributions follows from the convergence 
of the covariance, deduced from (|124l) and (I125p . □ 
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A Existence and uniqueness 

In this section we sketch the proof of existence, uniqueness and continuity of 
the solution of the system ([5^ . 

Theorem 4. Let B be the Banach space of vectors {X, h) eRx _L-^(R) with the 
norm 



II {X,h) \\B:=^\X\^ + \\h\\l. 
Let us consider the following Cauchy problem with initial datum {Xo,ho) G B 



ft 

X{t)=Xo + b{t)+ I dsT{X{s),h{s)) 



h{t)^p]ho~l db(s)p'l_,ipx{s)- I dsT{X{s),h{s))pJ_^ipx(s), 



(127) 



where T.'B ^ K is bounded and globally Lipschitz; recall that ip is a probability 
density in the Schwartz class of test functions and (px — ^{x — X). 
Then for any (Xo^ha) G B there exists a unique solution to {12'1I^ ; such a 
solution, {X{t),h{t)), belongs to C(R+;-B) and is such that 

s\Yp E\\{X{t),h{t))\\\ < oo Vr>0. (128) 

te[o,T] 

Uniqueness holds in the following sense: if {X{t),h{t)) is another continuous 
solution satisfying il28\) . then 

pI sup \\{X{t),hit))-{X{t),hm\l=o] =1 VT>0. 
\te[o,T] / 

Proof. We prove existence by Picard iterations, uniqueness by the Gronwall 
Lemma and continuity by using Kolmogorov's criterion. For the time being p] 
is either (HH) or so 7 S (0, 1). 

Existence: construct the sequence {(X(")(t), /i("'(t))} such that {X^^\h["'') = 
{Xo,p]ho) and, for n > 1, 



X^"\t) = Xo + b{t) + / dsT(X("-i)(s),/i("-i)(s)) 
Jo 

h<^"\t)=ptho~ fdb(s)pUp>xi^-^Hs) - fdsp]_J{X(--'\s),h^--'\s)), 
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where we set h) :— T{X, h)(px', notice that for a suitable constant K > 1 
we have 

\TiX,h) \^+mX,h)\\l + yx\\l<K 
\T{X, h)-T{Y,g)\ + \\(3{X,h)-l3{Y,g)h + yx-VYh<mX, h)-{Y,g)\\B, 
for any {X, h) and {Y,g) in B. Hence 

moreover, by the Cauchy-Schwarz inequality, 



E 



< t / dsE 







for n>l. Similarly, 



-2tE / ds 



p7_JT(X(")(s),/i(")(s))-T(X("-i)(s),/i("-i)(s))] 



Being p] a probability density, and because |jp7<<5||2 < llpZllill'Plb, p] is contrac- 
tive on L^(M); therefore 



"'0 

Iterating we end up with 



which gives uniform convergence on compacts [0, T] of the sequence (X^"^ (t), /i^"^ (i)) 
to a limiting process, {X{t),h{t)). Such a process is therefore an J^j-adapted 
solution to (IT?7)) . 

Uniqueness: by what we have done so far, it is clear that one can find a suitable 
c{t) uniformly bounded on compacts such that if (^X{t), is another solution, 
then 

E\\ MO) - {xit)Mml < cit) f dsE\\ {x{t)Mt)) - {x(t)Mml, 

Jo 

hence uniqueness follows by the Gronwall Lemma; (|128p is then a consequence 
of continuity, which we are going to prove. 

Continuity: being b(t) a.s. continuous and /3{X^h) bounded, X{t) is a.s. con- 
tinuous. In order to prove continuity for h(t) we first need to prove that for any 

geL^m 

Iim||p75-ffll2 =0. 
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In fact, using the scaling property of the kernel and the Jensen inequality 
(weighted version), we get 

11^75-5112= dx dwpl{w){g{x-wt'^)- g{x)) 

Jr Ur 

<fdxfdn.pJM igix-.n-gix)f 
Jr Jr 

= / dwpJ{w)\\T^t''g - g\\l 
Jr 

where Tr, t G M, is the translation {TTg){x) = g{x — r). Let us study the 
integrand: 

\\Tr9 - 9\\i = cwfTg - gwi = [ da e-'«^ m - m) f 

Jr 

=> linit_>o \\Tu,t-fg — 3II2 — for ^-C- w and 

pj{w) \\T^t.g - g\\l < CpJ{w)\\g\\l e L\R), 

so wc can apply the dominated convergence theorem and conclude. 
We are left with the continuity of k{t) := h{t) — plho. 

-k{t + 6) + k{t) = I dh{s){p1^^_^-pj_,)^x{s) 
Jo 

pt+S 

dh{s) p1^s_^ipx{s) 

f dsT[X{s),h{s)){P: 
Jo 

rt+S 

+ dsr{X{s),h{s))p'J^s_^ipx{s)- 



+ 



+ 



■t+S-s - Pl-s)fx(s) 



From now on we treat the cases < 7 < ^ and ^ < 7 < 1 separately. 
Let us start with the superdiffusion: 



E II k{t + 5)- k{t) \\i< + A2 + ^3 + Ai), 



where 



:= E 



ds {pI+s-s- pI-s)'Px(s) 



A2 ~E 



i-t+S 



'Pt+S-s'fiXis) 



A. := E 



f 

Jo 



db{s) {p1+s-s - Pt-s)fx{s) 
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Aa := E 



t+5 



We need to estimate all the above terms: 



Ai < CE 
= CE 

^CE 

< CE 

< CE 



t n4 



dsl dxl / dzpj{z)[ip{x — z{s + S)'^) (p{x — zs"')] 



2\ 2 



-, 4 



ds I dzp{{z)\\ip^(^s+S)y - 'fzs-'h 
< Cf^S^^ 



ds / dzpj{z) \ z \ 6'^ 
having used the scaling property and (1551) 

A2<E' 



dxS / {p]+s-sVxis)) ds 



= ^ E\^J^ dsWp-yipxit+s-sMj <C5\ 

having used the Cauchy-Schwartz inequality and the contractivity. 
In order to find estimates on the last two terms, let us choose i^{x) — ^1+ 
so that V/ e L2(R) , |j/||4 < llV'-^llill/V'lll • Hence, via the Burkholder 
equality and again Cauchy-Schwartz, we get 

A3 < U-^\\lE f db{s)ij{p]^s_^-p]_J^xis) 
Jo 

<CE f ds[^{p1^,_^~pU)^x(s)\ 
Jo 

<Ctj\sEj^dx^{x + X{s)f [{p]+s-s-pts)^Ti^) 
<Ctil + E sup \X{u)\') f dsUipl^s-PlMl 

u6[0,t] Jo 
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having used tp{x + X)"* < (1+ | X p)-0^(a;). Let us look at the integrand: since 
tjj{x) < tpiy) + \/| x — y I, we have 



mpi+s-piMt<c\\ip:+s-p:)ii^^)\\ 



+C dx 



dyip^+si^ - y)- pHx - y))y/\ ^-y My) 



(129) 



(130) 



The first addend can be estimated similarly to what we have done for Ai , so we 
get 

\\{p:+s-pz){i^v)\\i<c5'''- 

for the second, after applying Cauchy-Schwartz on the integrand, we find 



<c jjxUj^ dy{pi^s - pl)i^ - y) 



I a; - y I 



dyip^.s-pDix-y)^ (y) 



<C(^J^ dzpjiz) \z\{{s + Sr - s-')^ ' II (pI+, - p:)^Y2 < C5'\ 



and we end up with 



A3<Ct^\l + E sup I X{u) p 5^''. 



For Ai, analogously. 



Ai<C5\l + E sup I X{u) Y / ds||?Ap>||* 

V tiG[0,T] / Jo 



<C5\1 + E sup I X{u) p 



X ^ ds|||V'¥'lll + ^ dx (^J^ dyp1{x-y)yj\x-y \ip{y) 
Now the integral on the second line is estimated from above by 

' d J 11^^11^ + f/ dzpz{z) \z\ \ iipViin, 







so that 



Ai<C5'^ [l + E sup I X{u) p 



Proving continuity in the subdiffusive case is slightly more delicate; let us write 



E\\k{t + 5)- k{tW2 < CiAi +A2 + A3 + Ai), 



38 



where N = N{'y) is to be specified in the following and 



For A2: 



Ai :=£ 



An := E 



A3 :=E 
A4 := E 



Jo 

rt+5 

I db{s) {pI+^_, - Pt-s)<Px{. 
Jo 



2N 



2N 



t+S 



db{s) p'l+s_,V>x{s) 



2N 



A2 < cs 



2N 



Jo 



2N 



so that we need N > j. 

For ^3: let us choose again tp{x) = \ x \ as an auxiliary function; then 

VA^ > 0, ll^-'ll^ < 00 and V/ € L^R) , < \\ip-^\\^ ■ Via 

N-1 N -1 

the Burkholder inequality, using ip'^'^{x + X) < C(l+ | X{u) \^)tlj^^{x) and 
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working as we did for ^3 we get 



^3 < CE 



f ds V'^ [{Pt+s-s - Pl-s)fX(s)] 

Jo 

KCt^'-'il + E sup \ X{u)n fds [ dx ^'''{x) lipids -PD^ri^) 
ue[o,t] Jo Jr 

<Ct^-\l + E sup \X{u)f) [ ds [ dx iPix) dTp';'ip 
ue[o,t] Jo Jr Js 



(131) 



< Ct^-\1 + E sup I X{u) 1^) 
ue[o,t] 



fdsf 
Jo Jr 



X 

.JV-l 



dx 



■27 



2JV 



Ci^-^(l + £ sup I X{u) r) 

■[tG[0,i] 



fdsj 

Jo Jr 



dx 



ip{x) 


ll 







fS+S 



du pIlp" 



Jo {s + d-uy-^y Jo (6--m)1-27 



du plip" 



2N 



<Ct^-\l + E sup I Xiu) 1^) 
ue[o,t] 



X < ds dx 



%p{x) I duplip" 



1 



1 



+ / ds I dx 



s+S 



ipix) I du 



(s + 5-u)l-27 (s_u)l-27 
2N 

plv" 



2N 



(s + (5-u)l-27 



<Ce'-'{l + E sup I |^^)[^3a+^3f,], 

ue[o,t] 



(132) 



where 



^3a = I ds I dx I du 
Jo Jr Jo 

ds dx Id 
Jo Jr Jo 



1 



1 



1 



2JV 



(s + 5-m)1-27 (s 
t r rs+5 



/Jr 



2N 



1 



2JV 



2JV 
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We claim that 



dx 



dx 



max _{plip"i:){x) 

{)<u<s+d 



2N 



< OO, 



max (p2(-)%A*</?")(2:) 

0<u<s+d 



2N 



< OO. 



Indeed, p'^ip"ip is continuous in u, so the maximum in (|133P i is attained at, 
say, u and ||p2¥'"'0I12JV — ^- The maximmn in (1133^ 9 is reached at the second 
extremum (s + 5), in fact 



dx 



dypZix - y)^/\ x-y\ ip"{y) 



2N 



<C dxi dzpjiz) \z\^ 



2N 



Therefore, 



A3a<C ds 



du 



with C{t) bounded on compacts and 



2N 



2'y)2N 



A3b<C I ds 



s+S 



du 



{s + S-uy-^-r 



2N 



In order to apply the Kolmogorov criterion we need 4Nj > 1 and (1 — 27)2iV > 1. 
For Ai and A4: 

2N 



Ai<C 



dsip{x){p. 



t+5-s 



pI-s)VX(s) 



2N 



< Ct'^-'E dx ds\ ^{x)ip]^s-s ~ pU)vx(s) r\ 

JM. Jo 

which is exactly (llSip . 



Ai<C{l + E sup \X{u)\^^)S^^-^ f 
tie[o,t] JQ 



2N^ 



with analogous calculations, the integrand on the right hand side is bounded, 
hence 

^4<C(l+£^SUp I |2^)j2iV. 

«e[o,t] 

To conclude, requiring 

1 



N > 



N>ii if 7 < i, 



continuity follows. 



□ 
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B Motivation 



In the introduction we have briefly discussed the choice of the operators of 
fractional differentiation and of the fractional Laplacian. In this Appendix, we 
want to show how the operators DJ and naturally arise in the context of 
anomalous diffusion and explain in some more detail the link with CTRWs. 
We want to determine an operator A s.t. 

dtpJ{x)^Ap'l{x) 
pliO) = So, 

with p^{t,x) enjoying the following three properties: 

[dxp]{x)^l, [dxpj{x)x^0 e ( dxp2{x) x^ ^ t^^ (133) 

(notice that for 7 = i we recover the diffusion equation with ^ = A). We 
recall that /, and / denote the Fourier, the Laplace and the Fourier-Laplace 
transform of the function /, respectively. 
By (|133p . the following must hold 

p^^(k) = l~^ct^^'k^ + o{k^) and 

p^p, k) = -- -^r(27 + 1) = -(1 - cip-'-'k'), 

p 2p^^+^ p 

where ci = ^cr(27+l). In definitions ([7]) and ([U the constant ci should appear; 
we just set it equal to 1 both for simplicity and not being interested, in this 
context, in estimating the "anomalous diffusion" constant. 
We can assume that the expression for fP{p, k) is valid in the regime p~'^'^k'^ « 
1. Actually, condition (I133p ^ is meant for an infinitely wide system and for long 
times. In other words, if A is the region where the particle moves, we claim that 

lim lim — ^* ^ = const. 

This means that we are interested in the case k « p. Of course one can in 
principle find an infinite number of functions s.t. p'^{p,k) = ^(1 — cie) for 
e = p^^^k'^. One possible choice is 

1 /i^ 1 

P^{Pik) = —— = /i-)-!— — — — = — ■ (134) 

p{\-\-Cie) p'^^ + [ciky p + cik^p^ 

which leads to an integro-differential equation and, when 7 = i, it coincides 

with the Fourier-Laplace transform of a Gaussian density. 

We now find the operator whose fundamental solution is p^{p, k). We have 

C{dtp^{-,k)){p) = -l + pp-'{p,k)^^cik^p'-''^p'<ip,k). 



42 



Let p = 27 — 1 and (f>p (t) = ; then we need to distinguish two cases in order 
to study the right hand side of the above equation: 
when < 7 < ^ one can easily check that 

which imphes that 

1 /•* „ p''{s,k) . 



If 

{LL,k)p}^'^'' is the Laplace transform of — / ds- 

r(27-l)7o ( 



(t- s)2-27' 

when ^ < 7 < 1, instead, a straightforward calculation shows that 

c[^Mp+^*p■^{k,■))]^p^^i,k)p-p 

so that 

is the Laplace transform of — — - — / ds ^ \ . 

r(27) dt Jq [t~ s)^ 

Finally, taking the inverse Fourier transform, we get that p'^{t,x) satisfies (O 
when < 7 < ^ and ([5]) when i < 7 < 1. Moreover, the explicit expression for 
p'Jix) holds true: by (jl34l) we get that 

Jr 2^ 

hence 

2vcr 

and now, by the inverse Laplace formula, we obtain (jl6p . Obviously, the ex- 
pression (IT6|) has been deduced after having chosen (|134p among all possible 
candidates for p'^ and this choice can now be justified in view of the link with 
CTRWs. 
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